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ABSTRACT 


A study on mass transfer aspects of a sieve tray bench scale extraction column has been 
conducted on acetone/toluene/water system. From experiments, stagewise acetone 
concentration profile has been obtained for both dispersed and continuous phase which 
has then been compared with nonequilibrium and equilibrium based simulations of 
CHEMSEP and equilibrium based simulations of ASPEN PLUS software. Both UNIFAC 
and UNIQUAC method for LLE have been used in the simulations. Sum of error square 
analysis gives excellent matching experimental results with CHEMSEP nonequilibrium 
simulations without any stage efficiency; CHEMSEP equilibrium and ASPEN 
equilibrium results remain far off from the experimental results for both UNIFAC and 
UNIQUAC methods. 

Additionally photographic measurements gives drop size distribution and other 
hydrodynamic parameters which have been later used for calculating binary mass transfer 
coefficient of acetone in both dispersed and continuous phase. 



Chapter 1 
INTRODUCTION 

Liquid-liquid extraction or solvent extraction , one of the most important unit operations in 
chemical engineering is used for separating the components of a solution. It is one of the 
classical methods in separation technology and finds applications in the chemical and 
petroleum industry, hydrometallurgy, biotechnology, nuclear technology, waste 
management, and other areas [C. Tsouris et ah, 1994]. The two liquid systems are 
immiscible or partially miscible and are introduced into contacting equipment where a 
liquid dispersion with high enough interfacial area for mass transfer is created. A simple 
countercurrent liquid-liquid extraction process is shown in Figure 1.1. In the process, 
solvent extracts the solute from the feed solution. The solute rich phase leaving the column 
is called as ‘extract’ and the other phase depleted of solute leaving the column is called 
‘raffinate’. 

The liquid-liquid extraction operation consists of the following steps. 

a) intimate contacting of the solvent with the solution containing the 
component to be extracted so that the solute will be transferred from the 
solution to the solvent and 

b) separation of two immiscible phases. 

FEED EXTRACT 


SOLVENT 

Fig 1 . 1 simple liquid-liquid extraction column 
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1.1 Need for liquid-liquid extraction 

In any case, solvent extraction will be useful either where direct methods fail or where, 
despite apparent disadvantages, it nevertheless provides a less expensive overall process 
than a competitive direct or chemical method. The case for the use of liquid-liquid 
extraction will depend upon its either accomplishing a separation that can not be achieved 
by other operations such as distillation, evaporation and crystallization, or effecting the 
separation more economically. Liquid- liquid extraction is more economic alternative to 
other separation process and has found immense applications in the separation of ; 

a) solutions of components having low relative volatility, especially 
when vacuum distillation is expensive. 

b) Solutions of close boiling and azeotropic forming components. 

c) Dissolved solute when evaporation may be impractical. 

d) Solutions of heat -sensitive components such as antibiotics. 

\ e) Components of differing chemical type whose boiling points may 

overlap as in the case of petroleum hydrocarbons. 

But extraction has disadvantages also and these are listed below: 

a) regeneration of solvent is required to make the process feasible and 
environmentally acceptable 

b) Mutual solubility of the solvent and the feed causes contamination 
and necessitates further purification. 

c) Sometimes solvents are environmentally hazardous. So precaution 
has to be taken for using the system containing hazardous solvent. 

In liquid-liquid extraction process intimate contact of feed mixture with the solvent and the 
separation of the resultant mixture into two layers are essential. The liquid-liquid 
extraction equipment may be classified, according to the construction and the operational 

characteristics into two broad groups, stagewise contactors(mixer - settler) and differential 
contactors(continuous contactors). 

In differential contactors the composition of the phases change continuously. The flow 
may be parallel or countercurrent. They are more compact for a given throughput and 
usually require little ground area as they are normally in the form of vertical columns. But 
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the design and scale up of such units are, in general, complicated as the two phases never 
reach eqilibrium. 

1.2 Unagitated Sieve plate Column 

The perforated plate column consists of multiorifice plate with or without downcomers. 
These plates eliminate backmixing of the continuous phase totally. The light liquid is 
generally dispersed into drops and flows through the perforations and the continuous phase 
flows across the plate. The heavy liquid also could be dispersed by turning the column 
upside down. The dispersed drops coalesce and settle under or over the plate until 
sufficient hydraulic head is built up to force the liquid through the perforations. Sieve-plate 
columns as a liquid-liquid extraction equipment are widely used in petroleum refining 
operations [Muthuravichandran et ah, 1989] because of their uniqueness in providing 
repeated coalescence and redispersion of drops and offering advantages of cross-flow of 
the continuous phase. 


3 



Chapter 2 

LITERATURE REVIEW 


2.1 Mass transfer study in sieve tray LLX column on ternary system 

Few literatures are available where mass transfer studies has been done in bench scale 
sieve tray column. Many authors have developed extraction column design model and 
efficiency model for sieve tray extraction column where they used common ternary system 
as toluene /benzoic acid/water [Rocha and Fair, 1986], diethyl ether/acetic acid/water 
[Pyle, 1950], toluene/acetone/water [Rocha and Fair, 1986; Seibert and Fair, 1988; 
Seibert and Fair, 1993; Nanoti et. al., 1989], MIBK/acetic acid/water [Rocha and Fair, 
1986], 1-butanol/succinic acid/water [Seibert and Fair, 1988], 
Kerosen/carbontetrachloride/water [Wichterlova et. al., 1986]. Benzene/acetone/water, 
hexane/acetone/water and heptane/glycerine/water [ Muthuravichandran et. al., 1989]. 
But acetone /toluene/water system which has been chosen for our present study does not 
have many literatures where purely experimental work on liquid-liquid extraction has been 
done and detailed stagewise analysis of components were reported. [Rocha and Fair, 
1986] studied acetone/toluene/water system in a 10 cm (4 inch) diameter sieve tray liquid- 
liquid extraction column to develope mass transfer efficiency model. [Seibert and Fair, 
1988] was used this system in packed and spray column for hydrodynamic stydies. In this 
present work mass transfer study has been done for acetone/toluene/water system and 
sample analysis for both phases have been obtained by using Refractive Index method 
[Trevor et. al., 1996 ; Mario et. al., 1991] 

2.2 Models for binary mass transfer coefficient and diffusivity 

In liquid-liquid extraction,mass transfer occurs either from continuous phase to dispersed 
phase or vice versa. In dispersed phase, drop formation, drop fall/rise and drop coalescence 
are the three important phenomena in mass transfer. Binary diffusivities in the infinite 
dilution and that for concentrated multicomponent mixture are used to estimate the 
Maxwell-Stefan type binary mass transfer coefficient for both the phases. All the 

correlations for calculating mass transfer coefficient involve hydrodynamic parameters. In 
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literature several correlations are available for both diffusivity and mass transfer 
coefficient and these are discussed below. 


2.3 Binary diffusivity at infinite dilution 

Several correlations are available in the literature for estimating infinite dilution 
diffusivity. Some of those correlations are given below: 

[Wilke and Chang, 1955] 

( 2 . 1 ) 

where, :B°. is the Maxwel-Stefan type diffusion coefficient of species i present in 

infinitely low concentration in species j, T is temperature, is molar mass of the jth 

species, is molar volume of species i at its normal boiling point, is the association 

factor for the species j (2.26 for Water, 1.9 for methanol, 1.5 for ethanol, 1.0 for 
unassociated component) 

[Tyn and Calus, 1975] suggested the following correlation 

( 2 . 2 ) 

where P is the parachor 

[Hayduk and Minhas, 1982] reported that their correlation gave less error for diffusion of 
a solute in water. This correlation is given below 

=1.55*1 O'* (2-3) 

[Siddiqi and Lucas ,1986] correlation also had a noticeably lower average error for 
diffusing solute through water. 

:B° = 9.89 * 10"' jujO.907 j^-0.45p.-0.2652. (2.4) 

2.4 Binary Diffusivity in concentrated mixture 

Similarly several proposed correlations are available in the literature for estimating 
diffusivity:©,^. in concentrated mixtures. Most methods for predicting diffusivity in 

concentrated solutions attempt to combine the infinite dilution coefficients for and 

3©°. as a simple function of composition. 
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[Caldwell and Babb, 1956] proposed a simplest expression which is given as 

By =XjBy + x^B° ( 2 . 5 ) 

[Vignes, 1966] suggested that the composition dependence of By can be expressed by a 
relation of the form By = (By ) (B°. ) (2.6) 

But the Vignes equation is less successful for mixtures containing an associating 
component (e.g., an alcohol). Several modifications of the Vignes relation have also been 
proposed. 

[Leffler and Cullinan, 1970] for example included viscosity in the correlation as follows: 

(2.7) 

For specific system any one of the above models may give the good results. 

[Dullien and Asfour, 1985] found that the diffusivity for regular solutions may be well 

represented by By / = (B° / Y' (B° / ) ’'' (2. 8) 

[ Wesselingh and Krishna, 1990 ] proposed the following model for the binary Maxwel- 
Stefan diffusivity. By = (B° (B° (2.9) 

2.5 Binary mass transfer coefficient 

Using the above binary diffusivities, different correlations for binary mass transfer 
coefficients have been developed. It has been already mentioned that mass transfer takes 
place in three stages- drop formation, drop rise/fall, and drop coalescence. According to 
[Angelo and Lightfoot, 1968] approximately 90% of the mass transfer appears to be due 
to drop rise stage and essentially all the remainder to drop formation stage. [Philhofer, 
1984] assumed mass transfer due to fonnation to be about 5% of the total transfer. But 
[Rocha, 1986] based their model on Murphree stage efficiency concept which involved 
the use of overall mass transfer coefficients taking into account the mass transfer during 
drop fomation. 

2.5.1 [Handlos and Baron, 1957] proposed correlations for the mass transfer 
coefficients for transport within the both phases which are 

^|=0.00375C/,/(l + q,/M,) (2.10) 

A:^=0.725(Re)f^^(l-(/.,)C/,(S'c);^-^^ (2.11) 

where, is the slip velocity and is given by 
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+ ( 2 . 12 ) 
and M^are superficial velocity of dispersed and continuous phases, 0 ^ is 
fractional dispersed phase hold up. 

kfj , not being a function of implies that it is same for all binary pairs in the 
multicomponent system . 

2 . 5.2 [Skelland and Conger, 1973] have developed correlations for three different mass 
transfer stages in dispersed phase. During the drop rise (or fall) they have reported that 
selection of the appropriate correlations for the dispersed and continuous phase coefficients 
require knowledge of whether drops of the relavent size are internally stagnant, circulating, 
or oscillating. The correlations are shown below: 

2.5.2.1 Drop formation 

k‘ =0.0423(</„ lt,)(Ul ]■"■* (2-13) 

k; =0.386(S,;/(^)»"(p.r/Apj:l,M.)"’to//rf„)"'“ (2-14) 

where and are the densities of dispersed and continuous phase, Ap is density 
difference 

is the sautermean drop diameter, tj- is the drop formation time which is given 
by t,=(jtl6){dJI{inlA)dlV„} (2.15) 

Uff is the nozzle velocity, is the nozzle diameter and 7 is interfacial tension of 
the system. 

2.5.2.2 Drop coalescence 

k! = 0.173(J„, lyr^\Ult^ (2.16) 

kl = 5.959*10- (3©; It^ripJJl lgpy-^^\dlp,p,Ullpjr^^ (2.17) 

2.5.2.3 Drop fali/rise 

Internal circulation of drops in the commercial systems takes place with increasing the 
drop size. At lower size, drops are generally stagnant. Sustained oscillations of the drops 
apparently begin when the Weber number (We) reaches a value of 3.58 [Skelland and 
Conger, 1973]. 
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Stagnant drop 

kt /(0.5JJ] (2.18) 

where, is drop rise time, which can be expressed as ={Z^-h^)l 
= Compartment height, h, = mean height of the coalesced layer 
kl = (2.19) 

Circulating drop 

ki = 31.4(©' lyf" (2.20) 

=0.725(J„C/,p,/M,)^"(p./p,S,P^”!7,(l-,f,) (2.21) 

Oscillating drop 

= 0.32(©, ldJ{AB,.t^ ldjr\d^p,pjp,r\y^pl /gp^Apr^ (2.22) 

k;. = iB,j /dJ[50 + 0.00S5id^U,p, / p,y \pj p.B^jr (2.23) 


2.5.3 [Rocha and Fair, 1986] have presented correlations of the mass transfer coefficients 
for the drop formation zone and for drop rise/fall zone; the drop coalescence zone is held to 
contribute essentially nothing towards the total mass transfer in their rate equations. This 
correlations have been modified by Lao et. al. (1989) . correlations are given below: 

2.5.3. 1 Drop formation 


k^ = (-6.0 + O.OlWe + 6.5(r,^ lr))^B^.U, / d,^ 

(2.24) 

kl = (-6.0 + omwe + 6.5(7„^ lY))pp, Id,, 

(2.25) 


where. We = {d^Ulp^ lyg^) , rrp/=35 dynes/cm 

2.5.3.2 Drop fall/rise 

k^ = (0.70 + 0.03We){0.00315U^ /(I + p, / pj) (2.26) 

kl = (0J0 + 0.20We)0J25{d^U,pJ py-^\p, /(p^B‘)r°-^\l-<i>,)U, 

(2.27) 

2.5.4 [Seibert and Fair, 1993] developed correlations for drop formation and 
fall/rise,based on laboratory experiments on the system Toluene/AcetoneAVater. But as 
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like [Rocha and Fair, 1986] they also neglected the coalescence stage in their mass 
transfer calculations. 

2.5.4.1 Drop formation 


k^ =1.3(^)''' 

Ktf 

(2.28) 

CA 

II 

(2.29) 

2.5.4.2 Drop fall/rise 


if kl = o.mu,[\+Sc,r 

(2.30) 

<6, 0.00375(7, /(I + /ij ft,) 

(2.31) 

k;=0.69S(%Re°fScr(l-(j,,) 

^VS 

(2.32) 


[Debjit et al., 2000] developed a parallel mass transfer resistance model taking into 
accoxmt of all the three mass transfer stages which is as follows: 

(kDn. = C, {k ! ), + C, (k| ), + Q {k! ), (2.33) 


p refers to phase, d for dispersed phase and c for continuous phase 
where the correction factors can be expressed as 


r a a 

nr r' 


(2.34) 


where a’s are mass transfer area and C’s are area correction factors with respect to total 
area. 

2.6 Drop size measurement and distribution 

Drop size has a very important role on extraction efficiency. Extraction experiments 
[Tsouris et. al., 1994] revealed that mass transfer and drop behaviour are interrelated in 
liquid dispersions. In the case of solute transfer from the continuous phase to drops, 
smaller drop sizes are observed due to higher breakage rates. The smaller drops have 
longer residence times and, therefore, liie hold-up increases. In the other case of mass 
transfer, from the drops to the continuous phase, the coalescence rate is enhanced creating 
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larger drops and lower hold-up [Tsouris et al., 1994]. Mean drop size and drop size 
distribution measurements are done by photographic technique, followed by counting 
method. But the accuracy of the measurements depends on the choice of photographic site 
and proper optical adjustment. Drop may be non-spherical in shape. For non spherical 
drops, equivalent diameter can be estimated as follows [Vediyan et al., 1974] 

de = (dfd,r^ (2.35) 

where is the major axis and is the minor axis value. 

The mean size of the drops released at nozzle tips at low flow rates may be calculated by 
the correlation proposed by [Devotte and Chandrashekharan,1978]; 

d^(Y/Apgy'^ =23idlApglYr^^iUll2gd,r^^ (2.36) 

After counting number of drops by the photographic technique and their d^ values, the 
mean drop diameter which is called sauter-mean diameter ((i^^) and can be obtained as 
follows 

d.,=tdll±dl (2.37) 

1 1 

For the estimation of sauter-mean diameter they use 

d„ Hr /(Aat))" =1.59(17; ngd„)-*‘*' (2.38) 

where, d^ is the nozzle diameter , is the nozzle velocity 

According to [Laddha and Degaleesan, 1976] each compartment of a sieve plate column 
behaves like a spray column, and hence the mean drop size for the sieve plate column can 
be estimated by the equation number (2.38) 

From dimensional analysis, [Lao et. al., 1989] gave two correlations for the prediction of 
sauter-mean diameter. The correlations are given below; 


d.s 

/d^=l .591(fFe)'°“' 

for 

0< We <2 

(2.39a) 

d.s 

/d^=l .546(lFe)-®”' 

for 

2<We< 8.64 

(2.39b) 


The size of the drop increases with the increase in the nozzle diameter. To obtain the 
maximum drop diameter [Prabhu et al., 1967] recommended the following correlation. 

•i™. =1-38(7 /(Aa?)) (2.40) 
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The lower and upper limits of nozzle diameter which produced proper liquid jets with 
predictable distribution of drop sizes, were obtained by [Vediyan et al., 1972] and are 
given below. 

0.5(7 /(Apg)"' <d^< n(y /(Apg-)''^ (2.41) 

The lower limit is preferred for high interfacial tension systems and upper limit is for 
systems with low interfacial tension. 

[Nanoti et. al., 1984] used correlations for the prediction of drop diameter produced at the 
nozzle and they are given below. 

d^ld^= V' [2- 1 3(Ap / p , ) + exp(-0. 1 3Fr)] for E, < 0.4 (2.42a) 

djdf,= V'' [1 -24 + exp(-i^r )] for E, > 0.4 (2.42b) 

where E^ = Eotvos number 
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Chapter 3 


EXPERIMENTAL SETUP 

3.1 Design of Bench Scale Column 

3.1.1 Plate Design 

Choice of hole size 

Lower limit is preferred for high interfacial tension systems; that is for y greater than 20 
dynes/cm.The limits recommended by [Laddha et al.,1976] within which the hole size, 
dj^ may be chosen are given by 

0.5(y <d,<K{r /iApg)y'^ (3.1) 

Estimation of jet diameter, jet velocity and mean drop size 
[Treybal, 1963] suggests the following relationships to be used for obtaining the diameter 
of liquid jet dj issuing out of the perforation of the plate, depending upon the ratio of the 


hole diameter to the value of the drop characteristic group (y / Apgf'^ 

For <0.785 

d^/dj=l + 0.485[d^ KrApgy'^ f (3.2a) 

And for <7^ /(/ / ApgY^^ > 0.785 

d, Idj =0.12 + 1.51^ liY^PsY'"] (3.2b) 

Then nozzle velocity can be calculated from the following relation 

=1.69{djl d,f[{{y ! dj)}l{Q.5Up,+^Al\p^)f^ (3.3) 

It is also suggested by Treybal that the mean drop diameter of the dispersion produced by 
the jet may be taken as twice the jet diameter, i.e. 

(3.4) 

The actual sauter-mean drop size d„^ may be obtained from the nozzle velocity, , which 
[Vedaiyan, 1974 ] proposed 

d^liYl{Apg)r =\.592{Ull2gd^r-°^^^ (3.5) 
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Spacing of holes on the plate 

Since the drop size may be approximately 2 times the hole size at low flow rates, pitch 
should be at least 2 times of hole diameter. But 3 to 4 times of the hole diameter is 
recommended [Laddha, 1976]. In actual practice, the triangular pitch at distances which 
are multiples of hole size is often used. 

Dispersed phase head under the plate 

The total minimum liquid head h^ required for flow of the dispersed phase is given by 

h, = hfj + h^ + h^ (3.6) 

h, is also the thickness of the coalesced layer and is a measure of the static holdup, h^ is 

head required to overcome friction through perforations, h^ is head required to overcome 

interfacial effects and h^ is head required to overcome the effect of the flow of the 

continuous phase. According to [Bussalori et al.,1953] h^ can be expressed by 

K =(t/J[I-(S./S,)“lp,)/2CjAp (3.7) 

Sq is the area of all the perforations in the plate, S, is the total area of plate, is the 

orifice coefficient which can be calculated by [Laddha, 1976] 

Q=l-(0.71/logiV^J (3.8) 

where is the Reynolds number for the dispersed phase based on hole size and is 
given by =d^Uf,pJp^ (3.9) 

For estimation of h^ an empirical relationship is suggested by [Laddha, 1976] 

h, =(o.oirp;v“l/[Ap(</„)'‘] (3.10) 

where, / is in djmes/cm, ju^and are in centipoises, Ap is in lb/ ft ^ , and d^ and h^ are 
in inches. 

[Trebal, 1963] proposed the following equation for calculating h^ 

^P=4.5(i7i^p,/2Ap) (3.11) 

^down is the downcomer velocity, which can be calculated from the following equation 

t/.p«.=l-088(7Apg/p,^)>^"(t/^/2^rf^)-°“‘« (3.12) 

System properties have been taken at average values for composition range which were 
used in the design calculations, these average properties are given in the Table 3.1 
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Table 3.1 Properties of the system 


Property 

Toluene/ Acetone(28.20%Acetone+71.80%Toluene) 

Water 

Density( p) 

844 (p,) 

998 (p,) 

Viscosity(cp), ju 

0.56(p,) 

I-Ua^c) 

Interfacial tension 
(dynes/cm),7 

25 


All the design values obtained using equation (3.1) to (3.12) are given in the Table 3.2. 


Table 3.2 Design values of the column 


Parameter 

Equation used 

Value 

Hole or nozzle diameter, 

(3.1) 

3 mm 

Characteristic drop size, 
(y/Apgr^ 


4mm 

Jet diameter, dj 

(3.2a) 

2.36 mm 

Nozzle velocity, 

(3.3) 

13 cm/sec 

Mean drop 
diameter, (J^)„,„„ 

(3.4) 

4.72 mm 

Sauter-mean diameter, 

(3.5) 

6.63 mm 

Reynolds Number, 

^ P-d 

(3.9) 

814.8 

Orifice Coefficient, Cfj 

(3.8) 

0.756 


(3.7) 

16.50 mm 

\ 

(3.10) 

15.81 mm 

K 

(3.11) 

7.88 mm 

\ , coalesced layer thickness or 
static holdup 

(3.6) 

40.19 mm 


Plate spacing 

Normally in large scale unit 10 to 12 times the coalesced layer height is adopted mainly 
to provide flexibility for manual cleaning [Laddha,1976]. In actual practice for bench scale 
column plate spacing between 15 cm to 60 cm is employed. To find the sensitivity of plate 
spacing AcetoneAVater/Toluene system has been simulated in CHEMSEP software, 
developed by Horry Kooijman, Amo Haket and Ross Taylor, was used for 5 different 
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Acetone concentration in feed with varying tray spacing from 0.1m to 0.25m (0.1, 0.15, 
0.20 and 0.25). The design values of parameters that have been calculated from all the 
above equations (3.1) to (3.12), were used as input in this software. Desired flow rates for 
both phases were also used as input. Then % recovery of Acetone has been plotted against 
total number of stages (1,3,5,7,9,11,13 and 15). It has been observed from Fig A.1, Fig 
A.2, Fig A.3, Fig A.4, FigA.5 in Appendix that with increasing tray spacing recovery 
increases . But in CHEMSEP software sensitivity analysis of tray spacing has been done 
where sensitivity can be expressed as the following expression. 

^iray “ ^ ln(Re cov ery) / d In(trny) (3.13) 

The sensitivity study has been done for obtaining maximum acetone recovery in this 
extraction process. The sensitivity analysis plots Fig 3.1 and Fig 3.2 shows that after tray 
spacing of 0.2m there is no significant change in recovery . So tray spacing has been taken 
as 0.2m . 


3.1.2 Number of stages 

To calculate the height of the column it is very important to know the required number of 
stages. To find that, CHEMSEP software was used. The nonequilibrium column model has 

been chosen for running the Acetone/TolueneAVater system. The design values of 
parameters that have been calculated from all the above equations (3.1) to (3.12 
) and desired flow rates for both phases are used as input in this software. Desired flow 
rates for both phases were also used as input. Then different runs were carried out by 
varying stages from 1 to 15 (1,3,5,7,9,11,13 and 15) and varying solvent to feed molar 
ratio ranging from 3 to 7 (3,4,5,6 and 7). Different % Recovery of Acetone have been 
obtained for different run and the plots are shown in the Appendix. From these three plots, 
it has been observed that with increasing stages, recovery also increases. More than 80% 
recovery is obtainable using only 7 stages, then sensitivity analysis was studied where 
sensitivity is defined as 

= d ln(Re cov ery) I d hi{stage) (3.14) 

which are shown in the Fig 3.3 which indicates that after 8 stages there is no significant 
gain in recovery. So 8 stages have been taken as optimum number of stages. Standard 
column diameter has been taken as 10.2 cm [Rocha et al., 1978 & Seibert and Fair, 1993] . 
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A typical range of values of sieve tray geometries [Seibert and Fair, 1993 & 
Muthuravichandran et al., 1989] and design values of our column are given in the 
following Table 3.3 . 


Table 3.3 Sieve tray geometries of the column 


Parameter 

Typical Range 

This Column 

Hole diameter(cm) 

0.1-0.64 

0.102 

Fractional hole area 

0.03-0.08 

0.08 

Fractional downcomer area 

0.03-0.1 

0.038 

Pitch/hole diameter ratio 

3-4 

3 

Tray spacing (cm) 

15-64 

20 


3.2 Experimental Apparatus 

The chemical system chosen for the present study is Toluene/AcetoneAV ater. The liquid- 
liquid extraction was studied in a 10.2 cm diameter sieve tray column. The colmnn was 
made of three pieces of borosilicates glass section. They are 30cm, 34 cm and 98 cm in 
length. The longest section is fitted in the middle. Rubber packing separated the glass 
sections from each other. Aluminium flanges were fitted to the end of each glass section 
and the adjacent flange were connected by bolts. The column contained a central shaft 
fixed at the top and bottom. The peripheries of the downcomers were sealed by the M-seal. 
The shaft carried 7 number of sieve trays which were sealed by M-Seal with the glass wall 
to make leak proof. A cylindrical downcomer made of aluminium, was fixed near the 
periphery of each tray. A spiral shaped distributor was used to disperse the organic phase. 
To avoid corrosion the pipelines were made of S.Steel and Teflon, toluene and acetone 
used in this experiment were AR grade(99.50%), supplied by Thomas Baker chemicals 
Ltd., Kanpur. Twelve numbers of sample ports (stopcock made of glass) were fitted with 
the glass column. Six ports are rmder the tray and six are at the bottom of the downcomer 
as shown in the figure 3.4. Further details of the column and tray are given in Tables 3.4 


Table 3.4. Details of Plate and Downcomer 


Sr. No 

Item Description 

Value 

1 

Diameter of the column 

10.2 cm 
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2 

Effective height of the column 

178 cm 

3 

Diameter of the central shaft 

1.35 

4 

Thickness of the glass wall 

3 mm 

5 

Number of trays 

7 

6 

Plate spacing 

20 cm 

7 

Material of construction of plate 

Aluminium 

8 

Diameter of perforations 

3 mm 

9 

Total number of perforations 

100 

10 

Fractional free hole area of the plate 

0.08 

11 

Thickness of the plate 

1.5 mm 

12 

Diameter of the downcomer 

2.0 cm 

13 

Length of the downcomer 

12.5 cm 

14 

Number of sample port 

12 


3.3 Operating procedure 

The experimental set up is shown in figure 3.1 .Toluene and acetone mixture(dispersed 
phase) is pumped firom the storage tank ST2 to the bottom of the column through the 
distributor D. Water(continuous phase) is pumped to the top of the column. The pump 
used for continuous phase was a centrifugal pump(CP2) and a fire proof Pump(CPl) was 
used for dispersed phase. The distributor produced drops that moved countercurrently to 
the continuous phase. The flow rates of the both continuous and dispersed phase were 
regulated by needle valves(NVl and NV2). For accurate controlling of the flow rates, 
globe valves were used in the recycle line of both phases. Drops moved through the 
perforations of the trays and the continuous phase(water) moved through the downcomer 
attached with each trays. Total of 12 sample ports were located in the column for 
collecting samples for both phase and each stage. Dispersed phase samples were taken 
from the coalesced layer and continuous phase samples, firom the downcomer outlet. 

At the start of each run, column was filled with water approximately upto the continuous 
phase inlet. Then dispersed phase line was opened. The column was operated until steady 
state conditions were reached which was observed by the constant height of coalesced 
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layer under each trays. The column was operated for 45 minutes for each run.Usually 
steady state was attained after 20 minutes. After attaining the steady state of the column, 
all samples were collected (12 from sample ports , extract and raffinate) and were kept in 
airtight bottles. Then all samples were analysed in the following manner. 

3.4 Composition Analysis by Refractive Index 

For ternary composition analysis attempt has been done for Gaschromatograph. But due to 
the unavailability of this instrument binary analysis has been done by Refractometer. As 
solubility of toluene in water is very less, water presence in dispersed phase and toluene 
presence in continuous phase has been neglected. For calibration purposes different known 
concentrations of toluene /acetone mixture and acetone/water mixture were prepared. Then 
the refractive index of each sample was measured at a constant temperature (maintained by 
flowing water). The high precisioned refractometer used for the analysis was made by 
BAUSCH & LOMB company. The position of coincidence of crosswire with the interface 
gives the refractive index of the mixture. Using the refractive index value a standard plot 
was obtained of RI vs %concentration of acetone in the binary mixture which is linear. 
These plots have been shown in Fig 3. Sand Fig 3.9 and the calibration results are given in 
the Table 3.5 and Table 3.6. For each run, refractive index was measured of samples whish 
were previously collected. From the standard plots % acetone was measured for each 
sample in binary mixture. So for continuous phase % acetone was measured in 
water/acetone mixture and for dispersed phase % acetone were measured in 
toluene/acetone mixture. Calibration charts for Refractive Index of acetone/toluene mixture 
and acetone/water mixture are given in Table 3.5 and Table 3.6. 


Table 3.5 Calibration chart for acetone/water mixture 


% volume of 

Acetone volume 

Water volume 

Refractive Index 

Acetone in water 

(ml) 

» 

. . .(RI) 

10 

0.5 


1.33805 

20 

1.0 

4.0 

1.34408 

30 

1.5 

3.5 

1.34740 

40 

2.0 

3.0 

1.35380 

50 

2.5 

2.5 

1.35730 

60 

3.0 

2.0 

1.36020 

70 

3.5 

1.5 

1.36370 

80 

4.0 

1.0 

1.36860 
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90 

4.5 

0.5 

1.37220 


Table 3.6 Calibration chart for acetone/toluene mixture 


% volume of 

Acetone volume 

Toluene volume 

Refractive Index 

Acetone in Toluene 

(ml) 

(ml) 

(RI) 

10 

0.5 

4.5 

1.48778 

20 

1.0 

4.0 

1.47574 

30 

1.5 

3.5 

1.46060 

40 

2.0 

3.0 

1.44763 

50 

2.5 

2.5 

1.43658 

60 

3.0 

2.0 

1.42460 

70 

3.5 

1.5 

1.41850 

o 

oo 

4.0 

1.0 

1.39200 

90 

4.5 

0.5 

1.37910 


3.5 Drop size measurement 

The photographic method was used to study drop size and their distribution. For this 
purpose a rectangular view box (made of Perspex ,10cm x 10cm 15cm) was fitted at the 
cleanest part of the column through which photographs have been taken. 

Photographs were taken by using of a Pentax klOOO camera when column was operated 
under steady state condition. Especially this camera was used, as it is good for close 
photograph. For all sets of experiments lens used having power of 1.2 . The backside of the 
column was covered by black sheet to avoid any optical reflection. Three photographs 
were taken for each run. For measurement of drop sizes column shaft was used as 
reference scale. Minor and major axis were measured manually for non-spherical drops. 
The number of drops for the experiments of different concentration was counted . From the 

data of and , equivalent diameter was calculated from the equation number (2.35) 

and then sauter mean diameter was found out using the equation (2.37). 
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STl -Water storage tank, ST2-Feed storage tank, NV-Needle valve, 
GV-Gate valve,CPl-fire proof pump, CP2-centriiugal pump, R- 
Rotameter, Bl-Ball valve, ST3-Extract storage tank, ST4-Raffinate 
storage tank 


















Chapter 4 

RESULTS AND DISCUSSIONS 

4.1 Stagewise Concentration Distribution in the LLX Colunm 

The present work has been done on a 10.2 cm diameter bench scale LLX sieve tray 
column, which consists of 8 compartments. In this extraction operation, feed is 
acetone/toluene mixture and water has been used as solvent to extract acetone from the 
feed. In this operation, the solute acetone gets transferred from dispersed phase to 
continuous phase. All the 9 runs are isothermal and are combinations of 3 different solvent 
/ feed ratio and 3 different acetone concentration in the feed. After each steady state was 
achieved for each run, samples were collected from each stage, both for continuous and 
dispersed phase and then analysed by using a high precisioned refractometer. Hence 
composition of acetone in both phases have been measured. Thus a stagewise composition 
profile has been obtained for each run and plotted against the number of stage. 

All these results have been then compared with the simulator results. For this purpose two 
well known simulators CHEMSEP and ASPEN PLUS have been used to simulate the 
same runs on the same system. CHEMSEP has both equilibrium and non-equilibrium 
extractor model whereas Aspen Plus is a perfectly equilibrium model. In CHEMSEP non- 
equilibrium model, one can specify column specification and internal plate geometries . 
Using the experimental mole fraction and simulator mole fraction of acetone, residual sum 
of error square also has been calculated. Residual sum of error square is defined as 

For dispersed phase LSE’‘ = J - xf" 

1=1 

And, For continuous phase LSE'^ = L “ y predicted -y. 

1=1 

where N is the number of experimental data points and i is the component number. 

4.1.1 Experimental run specifications 

All the operation conditions of runs in the bench scale LLX column are given in the 
following Table 4.1. 


(4.1) 

(4.2) 
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Table 4.1 Flow rates and feed concentrations for experimental runs 


Run No. 

Feed flow 

Solvent flow 

S/F 

% Acetone in 


rate(gmol/sec) 

rate(gmole/sec) 


feed 

1 

0.085 

0.255 

3.0 

14.87 

2 

0.085 

0.340 

4.0 

14.87 

3 

0.085 

0.425 

5.0 

14.87 

4 

0.085 ' 

0.255 

3.0 

28.20 

5 

0.085 

0.340 

4.0 

28.20 

6 

0.085 

0.425 

5.0 

28.20 

7 

0.085 

0.255 

3.0 

40.25 

8 

0.085 

0.340 

4.0 

40.25 

9 

0.085 

0.425 

5.0 

40.25 


4.2 Comparison of experimental results with simulator results 

The experimental results obtained from LLX bench scale column has been compared with 
two well known software CHEMSEP and ASPEN PLUS. CHEMSEP software has both 
equilibrium and nonequilibrium extractor model, but ASPEN PLUS is an equilibrium 
model. Hence three models have been used and generated stagewise profiles have been 
compared with those of the experimental results. The runs given both of experiment and 
simulator have the same column, system and operating specifications already being shown 
in Table 4.1. 

Both UNIQUAC and UNIFAC property methods have been used to estimate the activity 
coefficient of the components in the acetone/toluene/water system. The method of 
estimation of activity coefficient is given in the Appendix in details. For using UNIQUAC 
method binary interaction parameters have been calculated in ASPEN PLUS and this 
values have been used in CHEMSEP software for simulation. For UNIFAC method this 
parameters are being estimated in the code inside. The results of comparison plots are 
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given in Fig 4.1(a) to Fig 4.9(b) for both UNIQUAC and UNIFAC property method. These 
overall plot shows a marked difference in the profile of equilibrium stage model of 
CHEMSEP and ASPEN PLUS With experimental results. The CHEMSEP equilibrium 
model has been run successfully for 6 cases out of 9 runs for UNIQUAC method and in 
UNIFAC method only in 4 runs code has been converged. The plots are showing that in 
maximum cases the experimental results are closely matching with CHEMSEP 
nonequilibrium results. 

With increasing feed concentration at constant solvent to feed ratio ASPEN equilibrium 
profile for UNIFAC is shifted downwards, at 28.20% acetone concentration in feed, it 
crosses the experimental data points as well as CHEMSEP nonequilibrium profile. But 
ASPEN profile for UNIQUAC does not show much of this effect. At lower feed 
concentration the CHEMSEP nonequilibrium profile for both UNIFAC and UNIQUAC 
overlap each other, it starts to move away at intermediate concentration and at higher feed 
concentration it produces a significant gap between them and in all cases the UNIFAC 
CHEMSEP nonequilibrium profile gives less error with experimental results compared 
with UNIQUAC CHEMSEP nonequilibrium results. 

At a fixed feed concentration with increasing solvent to feed flow rate , all profiles shifted 
downwards with respect to the experimental points. Also from the plots it has been 
observed that increase of solvent flow rate has a significant effect on mass transfer. With 
increase of solvent rate, the solute transfer zone in the column shifted from upper stages to 
lower stages when stage number is counted from the top of the column. The fact means 
that at lower solvent flow rate about all solute transfer occurs in some stages at the upper 
region of the column and this transfer shifted to the bottom zone in the column with 
increasing solvent flow rate. 

Both UNIFAC and UNIQUAC method have been used in all this simulation. But 
according to the behavior of the profiles, ASPEN profiles gives better sensitivity with 
changing the parameters (solvent to feed ratio and feed concentration) than CHEMSEP 
profiles. 

Sum of error square has been calculated for both CHEMSEP and ASPEN results for 
comparison with experimental results. These values are given in the Table 4.2 for 
dispersed phase and Table 4.3 for continuous phase. These value shows that for UNIFAC 
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method simulator results gives less error than UNIQUAC property method. Also the 
overall results gives the idea that CHEMSEP nonequilibrium model gives less sum of 
error square compared with other models. CHEMSEP equlibrium model gives m aximu m 
error for both in dispersed and continuous phase. 

Table 4.2 Comparison of error criteria for dispersed phase 


Run 

No. 

xlO'^ 

S! 

xlQ-' 

^ 3 ^ 

xlO'" 

5/ 

xlO"' 

S{ 

xlO"^ 

^ 3 ^ 

xlO”' 

1 

0.03 

roaigi 

15.2 

0.055 

13,24 

14.93 

2 



25.73 

0.136 

12.45 

28.66 

3 

0.282 

23 

30.8 

0.236 

7.21 

32.80 

4 

0.356 

6.26 

No convergence 

1.187 

14.31 

No convergence 

5 

2.83 

54.2 

No convergence 

0.236 

8.054 

222 

6 

1.03 

13 

136.4 

0.993 

1.450 

190.5 

7 

5.65 

31 

No convergence 

1.446 

7.780 

No convergence 

8 

4.91 

103.83 

No convergence 

0.603 

14.43 

No convergence 

9 

4.55 

226 

No convergence 

1.510 

32.119 

555 

Total 

19.766 

465.14 


4.979 

111.043 



Table 4.3 Comparison of error criteria for continuous phase 


Run 

Si 

S! 

S! 

■ 5 / 

^3^ 

s{ 

No. 

xlO"' 

xlQ-' 

xlO'^ 

xlO'" 

xlO'' 

xlO-' 

1 

0.002 


0.412 

0.003 



2 

0.001 


0.693 




3 

0.011 

0.161 


0.011 


32.80 

4 

0.082 

0.860 

No convergence 

■eeeshi 


No convergence 

5 

0.421 

0.308 

No convergence 

■eesshi 


5.36 

6 

0.026 

0.765 

1.266 

■EE&B 

0.115 

2.387 

7 

0.14 

2.61 

No convergence 


0.777 

No convergence 

8 


0.748 

No convergence 

■eeesh 

msmm 

No convergence 

9 


1.96 

No convergence 



5.271 

Total 





3.913 



Where, 
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cq — 

Oj chem(nonequilm) ^ ^2 Aspeniequilm)^ ^3 chem{equilm) ^ 

S^ — 5 ^2 “ Aspen(equilm) ^3 “" chem{equilm) 

This all are sum of error square for dispersed phase. For continuous phase this has been 
written as 

c9 _ A cpc — r c9 - r CF'" 

Oj chem(nonequilm) > ^2 Aspen (equilm) > ^3 chemiequilm) > 

and 

4.3 Evaluation of mass transfer coefficients 

4.3.1 Estimation of the Drop Size Distribution 

3 sample photographs that have been taken in the present study during the 9 runs are given 
in Fig 4.10 Fig 4.1 1 and 4.12. These photographs of drop swarais were obtained when the 
column was operating under steady state conditions. Dropsize measurements were made 
directly from the processed negatives. Three exposures were made, depending upon the 
drop population in the column to obtain a representation sample of the drop swarms. From 
the processed negatives, horizontal and vertical diameter ( and d^) of drops were 
measured and to these values, the proportionate factor, was applied in order to translate the 
photographic image to true image. For different types of drops, equivalent diameter was 
obtained by using d^ and d^ values in the equation (2.35) and these are given in the Table 
4.4, Table 4.5, Table4.6, Table 4.7 and Table 4.8.This d, and d^ values for different runs 
are given in the Appendix. The number of drops varies in our present study from 70 to 118 
with 10% error. The Fig 4.13 shows that with increasing acetone concentration in feed the 
number of drops in one compartment in the column increased. From the equivalent 
diameter of each type drops, the sauter-mean drop size were calculated by using the 
equation (2.37) . 

For each run drops were counted by using the photographs taken in a fixed compartment 
of the liquid-liquid extraction column and based on the equivalent diameters of each type 
drops, drop size distribution plots have been obtained which are shown in the Fig 4.14, Fig 
4.15 , Fig 4.16, Fig 4.17, Fig 4.18 . 

From Fig 4.14, Fig 4.15 and Fig 4.16, where plots obtained for three different Acetone 
concentration in feed (14.87%, 28.20% and 40.25%) and for same solvent to feed ratio( 
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3.0 ). For 14.87% acetone feed with increasing solvent rate number of drops has been 
reduced but for higher concentration as for 40.25% acetone in feed, number of drops 
increases with inceasing flow rate. The plots also shows that with increasing feed 
concentration, the unimodal curve is shifted to the left. The distribution is nearly normal or 
gaussian with minimum deviation at low solute concentration in feed, but with increasing 
solute concentration in feed, this deviation increases. This distributions are closely similar 
with the distributions obtained by [Vedaiyan et. al., 1972]. Also plot shows that peak 
diameter of the drop is becoming smaller with increasing feed concentration. It is also 
noticeable that at higher concentration larger drops are produced. In these plots it has been 
seen that at 40.25% acetone concentration maximum equivalent drop diameter produced 
10.944 mm where for 14 87% acetone concentration feed gives 8016 mm maximum drop 
size. Intermediate feed concentration ( 28.20% ) gives 9.28 mm maximum drop size. For 
the each distribution sauter-mean diameter has been calculated by using the equation 
(2.37) and values are given in the Table 4.10. 

Table 4.4 Drop size distribution for Run 1 


Equivalent drop diameter, (mm) 

Number of drops 

% contribution in total 

2.00 

2 

2.32 

4.51 

3 

3.16 

5.08 

9 

9.50 

5.82 

33 

34.74 

5.90 

15 

15.80 



6.752 

15 

15.80 

6.93 

10 

10.53 

7.92 

6 

6.32 

8.02 

1 

1.05 

10.01 

1 

1.05 

Total 

95 
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Table 4.5 Drop size distribution for Run 3 


Equivalent drop 

diameter, (mm) 

Number of drops 

% contribution in total 

4.80 

2 

2.86 

5.32 

3 

42.86 

5.47 

21 

30.00 

5.81 

4 

5.71 

6.35 

32 

45.71 

7.42 

3 

4.29 

8.45 

3 

4.29 

9.66 

1 

1.43 

10.98 

1 

1.43 

Total 

70 



Table 4.6 Drop size distribution for Run 4 


Equivalent drop 

diameter, (mm) 

Number of drops 

% contribution in total 

4.51 

7 

14.04 

5.81 

37 

32.46 

5.90 

16 

14.04 

6.30 

1 

0.88 

6.62 

31 

27.19 

7.11 

8 

7.02 

7.55 

2 

1.75 

8.38 

1 

0.880 

9.28 

2 

1.75 

Total 

105 



29 













Table 4.7 Drop size distribution for Run 7 


Equivalent drop 

diameter, (mm) 

Number of drops 

% contribution in total 

3.20 

4 

3.40 

3.71 

5 

4.24 

4.64 

25 

21.20 

5.82 

21 

17.80 

5.90 

14 

11.86 

7.62 

34 

28.80 

9.66 

7 

5.93 

9.98 

2 

1.70 

10.94 

6 

5.10 

Total 

118 



Table 4.8 Drop size distribution for Run 8 


Equivalent drop 

diameter, (mm) 

Number of drops 

% contribution in total 

4.19 

14 

12.30 

5.32 

11 

9.65 

5.47 

16 

14.04 

6.75 

34 

29.80 

7.19 

9 

8.00 

7.57 

11 

9.65 

7.62 

15 

13.16 

8.85 

1 

0.80 

10.12 

3 

2.63 

Total 

114 
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4 . 3.2 


Hydrodynamics Calculations 


Coalesced layer thickness has been measured for 5 runs ( runl, run3, run4, run? and 
run8) after the steady state achieved for each run and these are given in the Table 4.2.all 
this values has been measured by scale. 

Table 4.9 Coalesced layer measurement 


Stage 

Thickness of coalesced layer, A, , cm 


Rxinl 

Rim3 

Run4 

Run? 

Run8 

1 

2.4 

1 

1.5 

1 . 

1.1 

2 

2.5 

1.5 

1.7 

1.5 

1.5 

3 

3.0 

1.9 

2.1 

1.9 

1.6 

4 

3.0 

2.0 

2.0 

2.0 

1.6 

5 

2.8 

2.1 

1.6 

2.1 

1.5 

6 

1.5 

2.1 

1.7 

2.1 

1.2 

7 

1.2 

1.6 

1.1 

1.6 

1.0 


The fractional operational holdup can be defined as the ratio of the volume of the 
dispersed phase droplets to the effective volume of the column(Muthuravichandran et. al., 
1989). 

<j>j = (Total holdup- Sj. Y, K )/(total volume- (4.3) 

I 1 

where n is the number of plates, iSj-is the total area of cross section of the column, and 

tK=iK+K+K+K.+K+ ) (4-4) 

1 

where and are individual heights of the coalesced layers under plates 1, 

2, 3, 4 and 5 and so on and n is the total number of plates. 

Then sauter-mean diameter for each run has been calculated from the equation number 

(2.37) which are given in the Table 4.10 and using the sauter-mean diameter values, 
interfacial mass transfer area a has been calculated by the following equation 

a = 6^, Id, ( 4 . 5 ) 
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The vak perational holdup, Sauter-mean diameter and interfacial mass transfer area 
a are given in the Table 4.10. 

4.3.3 Calculation of binary mass transfer coefficient 

For the above 5 runs slip velocity has been calculated by using the equation number 
(2.12) in Chapter2. Slip velocity in this equation is a function of dispersed phase velocity 
, continuous phase velocity and fractional holdup <j)j . Then binary mass transfer 

coefficient of solute in dispersed and continuous phase has been calculated for acetone 
/toluene pair in dispersed phase using the [Handles and Baron, 1957] equation (2.10) and 
acetone/water pair for continuous phase using the equation (2.11). These values are shown 
in the Table 4.10. For equation (2.11) binary diffusivity value for acetone/water pair has 
been taken from the literature[Rocha and Fair, 1986]. In the dispersed phase, binary mass 
transfer coefficient of acetone/toluene pair has been calculated from the equation (2.10) 
and this equation has no diffusivity term. 


Table 4.10 Hydrodynamic parameters and binary mass transfer coefficient 


Run 

no. 

xlO'^ 

(mm) 

(2 5 

(m^ ! m^) 

(m/sec) 

xl0“^ 

(m/sec) 

xlO’^ 

(m/sec) 

xl0“^ 

NL 

k‘‘ 

^ace’-tol 

xlO"' 

(m/sec) 

ace-mter 

xlO'" 

(m/sec 

) 

1 

7.76 

6.447 

7.222 

1.12 

1.46 

16.00 

93.80 

3.976 

2.92 

3 

5.70 

6.607 

5.176 

1.12 

2.44 

22.20 

133.34 

6.588 

4.17 

4 

7.40 

6.416 

6.893 

1.05 

1.46 

15.70 

92.00 

3.900 

2.90 

7 

14.00 

7.625 

11.02 

1.00 

1.46 

8.61 

59.70 

2.140 

1.90 

8 

10.7 

6.983 

9.2 

1.00 

1.95 

11.30 

71.73 

2.810 

2.72 


4.4 Experimental error analysis 

From the previous all results it has been said that the experimental results are well 
matching with CHEMSEP nonequilibrium model. But some experimental data points 
produced significant deviation with CHEMSEP nonequilibrium results. This is because of 
experimental error. In the LLX column all sample ports for continuous phase has properly 
located. But the sample port in second and third stage in dispersed phase could not been 
properly located as there was a chance of breaking of glass column . Hence, during the 
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collection of sample some aqueous mixture also has been mixed with dispersed phase 
sample. So due to the refractive index error, it has been created erroneous composition 
analysis. 

Though the Refractometer which has been used for composition analysis was high 
precisioned because the accuracy is 0.00001, which gives 0.008% error in dispersed phase 
and 0.004% error in continuous phase in composition measurement. 
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mole fraction 



Figure 4.1(a) acetone concentration in dispersed phase for 
Runi ( 14.87% acetone in feed; S/F=3.0) 



Figure 4.1(b) Acetone cocentration in continuous phase for 
Run1 (14.87% acetone in feed; S/F=3.0) 
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Figure 4.4(a) Acetone concentration in dispersed phase for 
Run4 (28.20% acetone in feed; S/F=3.0) 



Figure 4.4(b) Acetone concentration in continuous phase for 
Run4 ( 28.20% acetone in feed; S/F=3.0) 
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Figure 4.6(a) Acetone concentration in dispersed phase for 
Run 6 ( 28.20% acetone in feed; S/F=5.0) 
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Figure 4.6(b) Acetone concentration in continuous phase for 
Run6 (28.20% acetone in feed; S/F=5.0) 
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Figure 4.7(a) Acetone concentration in dispersed phase for 
Run? ( 40.25% acetone in feed; S/F=3.0) 



Figure 4.7(b) Acetone concentration in continuous phase for 
Run? ( 40.255 acetone in feed; S/F=3.0) 
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4. 1 1 drop size distribution for 28.20% acetone in feed and S/F=3.0 



4.12 drop size distribution for 40.25% acetone in feed and S/F=3.0 
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Chapter 5 
CONCLUSION 


Liquid-liquid extraction experiments have been done on a sieve tray bench scale column 
with acetone/toluene/water system. Acetone is the solute which gets transferred fron 
toluene phase to aqueoue phase. A stagewise acetone concentration has been obtainec 
which has then been compared with CHEMSEPand ASPEN simulator results. Based or 
the results following conclusions have been drawn. 

a) Design has been done for a sieve tray bench scale LLX column by using standarc 
equations available in the literatures and CHEMSEP software. In CHEMSEI 
sensitivity analysis has been used for fixing an optimal tray spacing and numbei 
of stages. 

b) Experimental results match well with CHEMSEP nonequlibrium results. Bui 
ASPEN and CHEMSEP equilibrium model give significant error witt 
experimental results. However both for UNIFAC and UNIQUAC LLE methods. 
CHEMSEP nonequilibrium gives significantly less sum of error square, indicating 
that experimental results match well with nonequilibrium rate based model. 

c) Experimental results indicate that for both dispersed and continuous phase 
UNIFAC method gives better results than the UNIQUAC. 

d) Drop size distribution has been obtained for different runs in the experiment and 
these results are comparable to those in the literatures. 
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Recommendation for future work 

1 . Experiment can be done on the same system acetone/toluene/water by varying 
plate spacing, hole number on the plate, hole size and dispersed phase flow rates. 

2. Nonisothermal operation can be done on the extraction colxunn and the results can 
be compared with the simulator results. 

3. One can use different systems in the experimental work as mentioned in the 
Chapter 2. 

4. The result that obtained from the present study can be compared with simulation 
results using UNIQUAC method by incorporating experimental binary interaction 
parameters for the same system. 
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APPENDIX 


A.1 Models for calculating activity coefficient in multicomponent mixture 

A.l.l(a) UNIQUAC (Universal Quasi-chemical ) method 

Activity coefficients in liquid mixtures can be calculated from a model which expresses 
the excess Gibbs energy of the mixture as a function of the composition. This UNIQUAC 
method contains two adjustable parameters per binary. UNIQUAC in applicable to 
multicomponent mixtures of nonpolar and polar liquids(including those that participate in 
hydrozen bonding) as encountered in typical chemical and petrochemical processes. 


Excess Gibbs energy can be expressed as follows 

{combinatorial) + g ^ {residual) (A. 1 ) 

g^ {combinatorial) ! RT = x, ln((^), /x,) + X2 hi((^2 ! ^2)'^ + 92-^2 ^’'^(^2 ^^2) 

(A. 2 ) 

g'' {residual) / RT = 9,Xjln[0, (A. 3 ) 

T21 =expH(u2i r^i =exp{-[(M 2 i-Wu)/^ 2 "]} (A. 4 ) 

(j) is the average segment fraction: 

0 , = (XiTQ/(Xjrj +X2T2)and (/»2 = (x2T2)/(Xir, +X2T2) (A. 5 ) 

For multicomponent system the above equations will be 

g^ {combinatorial) / RT = J) x,. ln( 0 ,. / x,. ) + (z / 2 )^ g^.x,. ln( 0 ,. / ) (A.6) 

/ i 

g^ {residual) ! RT = (A- 7 ) 

i J 

where 

Xji = exp- {[(w ji -Uy)l RT] } (A.8) 

and the activity coefficient for component i becomes 

In 7,. = ln(^,. /x,.) + {zl 2 )qi ln(0,. /0, ) + /,■ - (0,- IXi)(^Xjlj)-qj ln(20 jT: p) + ~ 

J j J j 

(A. 9 ) 

where 

lj={z/ 2 ){rj-qj)-{rj-l) 


(A. 10 ) 


This and ip^ can be expressed by the following equations 
e,=q.N,/-Zq,Mj^q,x,/^q^x, 

J j 

/ S r.Nj = r,.x,. / Y. rjXj 

J J 

A.l.l(b) UNIFAC (UNIQUAC Functional Group) method: 


(A. 11) 
(A.12) 


In the UNIFAC method, the combinatorial part of the UNIQUAC activity coefficients, is 
used directly. Only pure component properties enter into this equation. Parameters r. and 

q. are calculated as the sum of the group volume and area parameters, and 

and =Yvi Qk (A. 13) 

k k 

where , always an integer, is the number of groups of type k in molecule i. Group 
parameters and are obtained from the van der Waals group volume and surface 

areas F,^^and given by [Bondi, 1968]: 

i?, = /15.17 and Q, = 4, /(2.5.10^) (A.14) 

15.17 and 2.5.10’ are the normalization factors. 

The residual part of the activity coefficient, can be written by the following equation 

Inrf =Zv'[lnr,-lnri] (A.15) 

k 

the group activity coefficient is found from an expression as follows 

InF^ =!2Jl-ln(E®»<PmJ-E(0m<Pfen/E0h<P„„) (A-16) 

m m n 

©„ is the area fraction of group m,and is group interaction parameter which are 
expressed as 

©m and (A.17) 

n 

(Pmn= exp- -U„JI Rr\ = exp- / T) (A. 1 8) 

where is the mole fraction of group m in the mixture and is a measure of the 

energy of interaction between groups m and n . 
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T able A. 1 .2(a) UNIQUAC binary interaction parameter estimated by ASPEN 
PLUS at 18°Cfor acetone/toluene/water system 


Component, i 

Component, j 


a .. 

Acetone 

Toluene 

148.213 

-78.455 

Acetone 

Water 

358.800 

-43.100 

Toluene 

Water 

1006.34 

373.22 


Table A. 1 .2(b) UNIFAC binary interaction parameter estimated by ASPEN PLUS 

at 18°C for acetone/toluene/water system 


Group ID i/j 

ACH 

ACCH, 

CH, 

CH,CO 

H^O 

ACH 

0 

167.00 

-11.12 

25.77 

903.80 

ACCH, 

-146.80 

0 

-69.70 

-52.10 

5695.00 

CH, 

61.13 

76.50 

0 

476.40 

1318.00 

CHfiO 

140.10 

365.80 

26.76 

0 

472.50 

H^O 

362.30 

377.60 

300.00 

-195.40 

0 
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% Recovery of Acetone 


1.2 Chemsep simulation of acetone/toluene/water syste 
acing and solvent to feed flow rate: 



Figure A.1 Acetone recovery for 14.87% acetone in feed w 
tray spacing 






Total no of stage 

Figure A.3 Acetone recoveiy for 40.25% acetone in feed with 
different tray spacing 












% Recovery of Acetone 



Figure A.7 Acetone recovery for 28.20% acetone in feed 
different S/F 



Fisure A. 8 Acetone recovery for 40.25% acetone in feed with di 
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Figure A. 9 calibration plot for acetone/toluene 
m ixture(Refractometer) 
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Figure A. 10 Calibration plot of acetone/water 
m ixtu re ( Ref racto m ete r) 




A. 1.3 Drop distribution and calculation of equivalent diameter for different Runs: 

Table A. 1 Equivalent drop diameter for Runl 


Major axis(mm) 

Minor axis(mm) 

Equivalent diameter(mm) 

2.0 

2.0 

2.00 

4.8 

4.0 

4.51 

6.4 

3.2 

5.08 

6.4 

4.8 

5.82 

8.0 

'3.4 

5.90 

8.0 

4.8 

6.752 

7.2 

6.4 

6.93 

8.8 

6.4 

7.92 

9.6 

5.6 

8.02 

11.2 

8.0 

10.01 


Table A.2 Equivalent drop diameter for Rim3 


Major axis(mm) 

Minor axis(mm) 

Equivalent diameter(mm) 

4.8 

4.8 

4.80 

5.6 

4.8 

5.32 

6.4 

4.0 

5.47 

6.4 

4.8 

5.81 

8.0 

4.0 

6.35 

8.0 ^ 

6.4 

7.42 

11.2 

4.8 

8.45 

11.2 

7.2 

9.66 

14.4 

6.4 

10.98 
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Table A.3 Equivalent drop diameter for Run4 


Major axis(mm) 

Minor axis(mm') 

Equivalent diameter(mm') 

4.8 

4.0 

4.51 

10.4 

4.0 

5.81 

8.0 

3.2 

5.90 

11.5 

4.8 

6.30 

7.2 

5.6 

6.62 

8.0 

5.6 

7.11 

10.4 

4.0 

7.55 

9.6 

3.2 

8.38 

11.2 

6.4 

9.28 


T able A.4 Equivalent drop diameter for Run? 


Major axis(mm) 

Minor axis(mm) 

Equivalent diameterfmm) 

3.2 

3.2 

3.20 

4.0 

3.2 

3.71 

5.6 

3.2 

4.64 

6.4 

4.8 

5.82 

8.0 

3.2 

5.90 

9.6 

4.8 

7.62 

11.2 

7.2 

9.66 

14.4 

4.8 

9.98 

12.8 

8.0 

10.94 
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Table A. 5 


Equivalent drop diameter for Run8 


Major axis(mm) 

Minor axisfmm') 

Equivalent diameter(rntn) 

4.8 

3.2 

4.19 

5.6 

4.8 

5.32 

6.4 

4.0 

5.47 

8.0 

4.8 

6.75 

8.8 

OO 

7.19 

8.8 

5.6 

_ ^ 

7.57 

9.6 


7.62 

10.40 

6.4 

8.85 

13.6 

5.6 

10.12 


A.1.4 Sample analysis by Refractometer: 

Table A.6 Dispersed phase sample analysis for Runl 


Stage 

RI 

% Acetone 

1 

1.49878 

4.48 

2 

1.49710 

6.35 

3 

1.49607 

7.49 

4 

1.49401 

9.73 

5 

1.49306 

10.76 

6 

1.49175 

12.16 

7 

1.49070 

13.26 


Table A.7 Continuous phase sample analysis for Runl 



RI 

% Acetone 

2 

1.3384 

0.89 

3 

1.33906 

1.41 

4 

1.33967 

1.87 

5 

1.34031 

2.37 

6 

1.34116 


7 

1.34209 

3.82 












Table A.8 


Dispersed phase sample analysis for Run! 


Stage 

RI ~ 

% Acetone 

1 

T 49933 

3.97 

2 

1.49775 ■ ■ 

5.80 

3 

1.49720 

6.41 

4 

1.49585 

7.97 

5 

1 .49495 ~ 

8.98 

6 

1.49329 

10.82 

7 

1.49099 

12.95 


Table A.9 Continuous phase sample analysis for Run2 


Stage 

RI 

% Acetone 

2 

1.33794 

0.59 

3 

1.33836 

0.89 

4 

1.33894 

1.32 

'5 

1.33957 

1.79 

6 

1.34023 

2.30 

7 

1.34086 

2.80 


Table A. 1 0 Dispersed phase sample analysis for Run3 


Stage 

RI 

Vo Acetone 

1 

1.50000 

3.13 

2 

1.49906 

4.28 

3 

1.49848 

4.96 

4 

1.49661 

7.10 

5 

1.49496 ! 

8.96 

6 

1.49291 

11.24 

7 

1.49216 

\2.01 


Table A. 1 1 Continuous phase sample analysis for Run3 


Stage 

RI 

Yo Acetone 

2 

1.33783 

0.52 

3 

1.33842 

0.93 

4 

1.33840 

0.92 

5 

1.33880 

1.22 

6 

1.33958 ' 

1.80 

7 

1.34025 

2.32 











Table A. 1 2 Dispersed phase sample analysis for Run4 


Stage 

RI 

% Acetone 

1 

1.49225 

11.97 

2 

1.49043 

13.95 

3 

1.48745 

17.13 

4 

1.48388 

20.83 

5 

1.48062 

24.13 

6 

1.47858 

26.14 


1.47802 

26.70 


Table A. 1 3 Continuous phase sample analysis for Run4 


Stage 

R[ 

% Acetone 

2 

1.33978 

1.95 

3 

1.33326 

2.87 

4 

1.34234 

4.03 

5 

1.34320 

4.78 

6 

1.34410 

5.60 

7 

1.34430 

5.78 


T able A. 1 4 Dispersed phase sample analysis for Run5 


Stage 

RI 

% Acetone 

1 

1.49189 

■■■■BSSSHHIHl 

2 

1.48856 


3 

1.48554 

18.58 

4 

1.48166 

22.44 

5 

1.48085 ^ 

23.24 

6 

1.47867 

25.34 

7 

1.47710 

26.80 


Table A. 1 5 Continuous phase sample analysis for Run5 


Stage 

RI 

% Acetone 

2 

1.33998 

2.11 

3 

1.34138 

3.22 

4 

1.34220 

3.92 

5 

1.34312 

4.17 

6 1 

1.34354 

5.08 

7 ' 


5.52 
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Table A. 1 6 Dispersed phase sample analysis for Run6 


Stage 

RI 

% Acetone 

i 

1.49430 

9.71 

2 

1.49290 

11.27 

3 

1.49026 

14.14 

4 

1.48635 

18.29 

5 

1.48308 

21.64 

6 

1.48068 

24.06 

? 

1.47930 

25.43 


Table A. 1 7 Continuous phase sample analysis for Run6 


Stage 

RI 

% Acetone 

2 

1.33825 

0.82 

3 

1.33965 

1.85 

4 

1.34039 

2.43 

5 

1.34121 


6 

1.34179 

3.57 

7 

1.34242 

4.10 


Table A. 1 8 Dispersed phase sample analysis for Run? 


Stage 

RI 

% Acetone 

1 

1.48970 

24.88 

2 

1.48070 

28.50 

3 

1.48340 

30.30 

4 

1.48700 

33.00 

5 

1.48733 


6 

1.48800 


7 

1.48890 

38.40 


Table A. 1 9 Continuous phase sample analysis for Run? 


Stage 

RI 

% Acetone 

2 

1.33770 

2.95 

3 

1.33750 

3.44 

4 

1.33733 

3.86 

5 

1.33725 1 

4.05 

6 

1.33720 ^ 

4.18 

7 

1.33710 

4.42 










Table A.20 Dispersed phase sample analysis for Run8 


Stage 

RI 

% Acetone 

1 

1.48804 

16.50 

2 

1.48155 

23.20 

3 

1.47924 

25.50 

4 

1.47310 

31.40 

5 

1.47139 

33.00 

6 

1.46701 

37.00 

7 

1.46612 

37.80 


Table A.2 1 Continuous phase sample analysis for Run8 


Stage 

RI 

% Acetone 

2 

1.34061 

2.60 

3 

1.34189 

3.65 

4 

1.34292 

4.53 

5 

1.34346 

5.00 

6 

1.34477 

6.22 

■ 7 

1.34527 

6.71 


Table A.22 Dispersed phase sample analysis for Run9 


Stage 

RI 

% Acetone 

1 

1.49001 


2 

1.48567 


3 

1.48223 


4 

1.47913 

25.60 

5 

1.47491 

29.70 

6 

1.47108 

33.30 

7 

1.46612 

37.80 


Table A.23 Continuous phase sample analysis for Run9 


Stage 

RI 

% Acetone 

2 

1.33966 

1.86 

3 

1.34044 

2.47 

4 

1.34173 

3.51 

5 

1.34285 ^ 

4.47 

6 

1.34371 

5.24 

7 

1.34433 

5.81 


&5 
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